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Abstract. If CDM particles decay and their lifetime is comparable to the age of the Universe, they can modify 
its equation of state. By comparing the results of numerical simulations with high redshift SN-Ia observations, we 
show that this hypothesis is consistent with present data. Fitting the simplest quintessence models with constant 
w q to the data leads to w q < —1. We show that a universe with a cosmological constant or quintessence matter 
with w q ~ — 1 and a decaying Dark Matter has an effective w q < — 1 and fits SN data better than stable CDM 
or quintessence models with w q > —1. 



There are at least two motivations for the existence of a Decaying Dark Matter (DDM). If R-parity in SUSY 
models is not strictly conserved, the LSP which is one of the best candidates of DM can decay to Standard 



Model particles Banks et al. 1995. Violation of this symmetry is one of the many ways for providing neutrinos 
with very small mass and large mixing angle. 

Another motivation is the search for sources of Ultra High Energy Cosmic Rays (UHECRs)(see Yoshida & 
Dai 1998 for review of their detection and Blandford 1999| and Bhattacharjee & Sigl 1998] respectively for 
conventional and exotic sources). In this case, DDM must be composed of ultra heavy particles with Mdm ~ 
10 22 — 10 25 eV\ In a recent work Ziaeepour 2000 we have shown that the lifetime of UHDM (Ultra Heavy Dark 
Matter) can be relatively short, i.e. r ~ 10 — IOOto where To is the age of the Universe. Here we compare 
the prediction of this simulation for the Cosmic Equation of State (CES) with the observation of high redshift 
SN-Ia. 



For details of the simulation we refer the reader to Ziaeepour 200C . In summary, the decay of UHDM is assumed 
to be like the hadronization of two gluon jets. The decay remnants interact with cosmic backgrounds, notably 
CMB, IR, and relic neutrinos, lose their energy and leave a high energy background of stable species e.i. e^, 
, u, P, and 7. We solve the Einstein-Boltzmann equations to determine the energy spectrum of remnants. 



Results of Ziaeepour 200C show that in a homogeneous universe, even the short lifetime mentioned above can 
not explain the observed flux of UHECRs. The clumping of DM in the Galactic Halo however limits the possible 
age/contribution. These parameters are degenerate and we can not separate them. For simplicity, we assume 
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Fig. 1. Energy density of the Universe. Solid line Vf£ = £!a = 0.7 and stable DM; dashed line the same 
cosmology with r = 5to; dash dot line A = and stable DM; dot line A = and r = 5tq. Dependence on the 
mass of DM is negligible. 



that CDM is entirely composed of DDM and limit the lifetime. Fig. [j] shows the evolution of energy density 
T 00 (z) = p(z) at low and medium redshifts in a flat universe with and without a cosmological constant. As 
expected, the effect of DDM is more significant in a matter dominated universe i.e. when A = 0. For a given 
cosmology, the lifetime of DDM is the only parameter that significantly affects the evolution of p. For the same 
lifetime, the difference between Mdm = 10 12 eU and Mdm = l0 24 eV cases is only w 0.4%. Consequently, in 
the following we neglect the effect of the DM mass. For the same cosmological model and initial conditions, if 
DM decays, matter density at z = is smaller than when it is stable because decay remnants remain highly 
relativistic even after losing part of their energy. Their density dilutes more rapidly with the expansion of the 
Universe than CDM and decreases the total matter density. Consequently, relative contribution of cosmological 
constant increases. This process mimics a quintessence model i.e. a changing cosmological constant Peebles & 
Ratra |l988[ Zlatev, Wang & Steinhardt |1999| (see Sahni & Starobinsky |1999| for recent review). However, the 
equation of state of this model has an exponent w q < — 1 which is in contrast with the prediction of scalar field 
models with positive potential (see appendix for an approximative analytical proof). 

The most direct way for determination of cosmological densities and equation of state is the observation of 
SN-la's as standard candles. It is based on the measurement of apparent magnitude of the maximum of SNs 



lightcurve Perlmutter et al. 1997 1999, Riess A. et al. 1998. After correction for various observational and intrinsic 
variations like K-correction, width-luminosity relation, reddening and Galactic extinction, it is assumed that 
their magnitude is universal. Therefore the difference in apparent magnitude is only related to difference in 
distance and consequently to cosmological parameters. 

The apparent magnitude of an object m(z) is related to its absolute magnitude M: 



m(z) = M + 25 + 51og£> L 



(1) 



where Dl is the Hubble-constant-free luminosity distance: 



dz' 
o W) 



sinh(a;) Qn> 0, 

<%) = x n R = o, 

sm(x) fin < 0. 



(2) 



(3) 



E(z) 



H{z) 



SttG, 



T ou (z) + 



(4) 
(5) 



Here we only consider flat cosmologies. 



We use the published results of the Supernova Cosmology Project, Perlmutter et al. 1999 for high rcdshift and 



Calan-Tololo sample, Hamuy et al. 1996 for low redshift supernovas and compare them with our simulation. 



From these data sets we eliminate 4 SNs with largest residue and stretch as explained in Perlmutter et al. 1999 
(i.e. we use objects used in their fit B). 

Minimum-x 2 fit method is applied to the data to extract the parameters of the cosmological models. In all fits 
described in this letter we consider M as a free parameter and minimize the \ 2 with respect to it. Its variation 



in our fits stays in the acceptable range of ±0.17, Perlmutter et al. 1997 



We have restricted our calculation to a range of parameters close to the best fit of Perlmutter et al. 1999 i.e. 
2.38 x 10~ n ^ pa = ^ 3.17 x 10~ 11 eV 4 . The reason why we use pa rather than Qa is that the latter 
quantity depends on the equation of state and the lifetime of the Dark Matter. The range of pa given here 
is equivalent to 0.6 ^ fl e A q ^ 0.8 for a stable CDM and Ho = 70 km Mpc -1 sec -1 (we use $l A q notation to 
distinguish between this quantity and real J7a). 

Fig.|| shows the residues of the best fit to DDM simulation. Although up to 1-cr uncertainty all models with 
stable or decaying DM with 5to < r < 50to and 0.68 < QT* < 0.72 are compatible with the data, a decaying 
DM with t ~ 5to systematically fits the data better than stable DM with the same A 9 . Models with A = are 
ruled out with more than 99% confidence level. 

In fitting the results of DM decay simulation to the data we have directly used the equation (||) without defining 
any analytical form for the evolution of T 00 (z). It is not usually the way data is fitted to cosmological models 
Perlmutter et al. 1997, 1998, Garnavich et al. 1994 Consequently, we have also fitted an analytical model to the 
simulation for z < 1 as it is the redshift range of the available data. It includes a stable DM and a quintessence 
matter. Its evolution equation is: 



8ttG 



F 2 (z) = ^ (T 00 + n ^ + 1) 3K +1)) _ 



(6) 



The term T s ° t ° is obtained from our simulation when DM is stable. In addition to CDM, it includes a small 
contribution from hot components i.e CMB and relic neutrinos. For a given Q A q and r, the quintessence term is 
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Fig. 2. Best fit residues with Vf£ = 0.7, r = 5tq. It leads to f^A = 0.73. The curves correspond to residue for 
stable DM with £f£ = fl\ = 0.7 (doted); A = and r = 5r (dashed); A = 0, stable DM (dash-dot). 



fitted to T 00 - T°° + The results of this fit are £l q and tu g which characterize an equivalent quintessence 
model for the corresponding DDM. The analytical model fits the simulation extremely good and the absolute 
value of relative residues is less than 0.2%. Results for models in the 1-cr distance of the best fit is summarized 
in Table 0. 

In the next step, we fit an analytical model to the SN-Ia data. Its evolution equation is the following: 

h 2 (z) = ^f((i - n q )( z + 1) 3 + n q (z + i)3(«.+i)). (7) 

The aim for this exercise is to compare DDM equivalent quintessence models with the data. 
Fig.|| shows the \ 2 of these fits as a function of w q for various values of £l q . The reason behind using \ 2 rather 
than confidence level is that it directly shows the goodness-of-fit. As with available data all relevant models are 
compatible up to 1-cr, the error analysis is less important than goodness-of-fit and its behavior in the parameter- 
space. 

Models presented in FigJ| have the same fl q as equivalent quintessence models obtained from DDM and listed 
in Tab.[j]. These latter models are shown too. In spite of statistical closeness of all fits, the systematic tendency 
of the minimum of x 2 to w q < — 1 when fl q < 0.75 is evident. The minimum of models with fl q > 0.75 has 
w q > —1, but the fit is worse than former cases. Between DDM models, one with fl q = 0.71 is very close to the 
best fit of (^) models with the same Q q . Regarding errors however, all these models, except Q q = 0.8 are 1-cr 
compatible with the data. 

One has to remark that Q q and w q are not completely independent (see Footnote [l]) and models with smaller 
1 The exact equivalent model should have the same form as ([?]). However, it is easy to verify that in this case, the 
minimization of \ 2 has a trivial solution with w q = — 1, Q q = 0. Only one equation remains for non-trivial solutions and 
it depends on both w q and Q, q . Consequently, there are infinite number of solutions. 

The model we have used here generates a very good equivalent model to DDM with less than 2% error (Because CDM 
and quintessence terms are not fitted together, Q is not exactly 1). 



Table 1. Cosmological parameters from simulation of a decaying DM and parameters of the equivalent 
quintessence model. H is in km Mpc -1 sec -1 . 
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Fig. 3. x 2 -fit of models defined in (|?|) as a function of w q for 9 = 0.67 (dashed), fl q — 0.69 (dash-dot), 
£l q = 0.71 (solid) and fl q = 0.8 (dotted). The x 2 of equivalent quintessence models to DDMs with r = 5ro and 
same fi g is also shown. Except Q g = 0.8 model, others are all the best fit to DDM. For Q q = 0.8, a stable DM 
fits the data better, but the fit is poorer than former models. 



£l q and smaller w q has even smaller \ 2 - ln fact the best fit corresponds to il q = 0.5, w q = —2.6 with x 2 — 61.33. 
The rejection of these models however is based on physical grounds. In fact, if the quintessence matter is a 
scalar field, to make such a model, not only its potential must be negative, but also its kinetic energy must be 
comparable to the absolute value of the potential and this is in contradiction with very slow variation of the 
field. In addition, these models are unstable against perturbations. It is however possible to make models with 



w q < — 1, but they need unconventional kinetic term Caldwell 1999 



These results are compatible with the analysis performed by Garnavich et al. 1998. However, based on null 
energy condition Wald 1984, they only consider models with w q ^ —1. This condition should be satisfied by 



non-interacting matter and by total energy-momentum tensor. As our example of a decaying matter shows, a 
component or an equivalent component of energy- momentum tensor can have w q < — 1 when interactions are 
present. 

In conclusion, we have shown that a flat cosmological model including a decaying dark matter with r ~ 5to 
and a cosmological constant or a quintessence matter with w q ~ —1 at z < 1 and Q q ~ 0.7 fits the SN-Ia data 
better than models with a stable DM or w q > — 1. 

The effect of a decaying dark matter on the Cosmic Equation of State (CES) is a distinctive signature that can 
hardly be mimicked by other phenomena, e.g. conventional sources of Cosmic Rays. It is an independent mean 
for verifying the hypothesis of a decaying UHDM. In fact if a decaying DM affects CES significantly, it must 
be very heavy. Our simulation of a decaying DM with M ~ l0 12 eV and r = 5to leads to an over-production 
by a few orders of magnitude of 7-ray background at E ~ 10 9 — 10 n eV with respect to EGRET observation 



Sreekumar et al. 1998 Consequently, such a DM must have a lifetime much longer than 5tq. However, in this 
case it can not leave a significant effect on CES. 

The only other alternative for making a quintessence term in CES with w q sightly smaller than — 1, is a scalar 
field with a negative potential. Nevertheless, as most of quintessence models originate from SUSY, the potential 
should be strictly positive. Even if a negative potential or unconventional models are not a prohibiting conditions, 
they rule out a large number of candidates. 

Present SN-Ia data is too scarce to distinguish with high precision between various models. However, our results 
are encouraging and give the hope that SN-Ia observations will help to better understand the nature of the 
Dark Matter in addition to cosmology of the Universe. 

Appendix: Here we use an approximative solution of (|J) to find an analytical expression for the equivalent 
quintessence model of a cosmology with DDM and a cosmological constant. With a good precision the total 
density of such models can be written as the following: 
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, T - t 
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(8) 



We assume a flat cosmology i.e. £Im + fiA = 1 (ignoring the hot part). p c is the present critical density. If DM 
is stable and we neglect the contribution of HDM, the expansion factor a(t) is: 
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Using (JsJ) as an approximation for when DM slowly decays, (|sj> takes the following form: 
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For a slowly decaying DDM, ar >> 1 and ( [l3| ) becomes: 
POO 



Pc 
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fi 9 (l + z) 3 ^ = fi A (l + ^-z(l + z) 3 lnC) 



Equation (16) is the definition of equivalent quintessence matter. After its linearization: 



i fi M (l + 4A)(l- V2A) 
»«=7,-l« I- 

It is easy to see that in this approximation w q < — 1 if f^A > 5. 



(13) 
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